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SUMMARY 

The vibrational behavior of cantilevered aircraft wings modeled as thin— walled beams and 
incorporating piezoelectric effects is investigated. Based on the converse piezoelectric effect, the 
system of piezoelectric actuators conveniently located on the wing yield the control of its 
associated vertical and lateral bending eigenirequenties. The possibility revealed by this study 
enabling one to increase adaptively the eigenfrequencies of thin-walled cantilevered beams could 
play a significant role in the control of the dynamic response and flutter of wing and rotor blade 
structures. 


INTRODUCTION 

The successful development of smart material systems technology [Rl] is likely to generate 
new avenues and concepts toward the design of the next generation of aeronautical and aerospace 
vehicles. In spite of the complexity and severity of environmental conditions to which these 
vehicles are likely to be exposed, they must be designed as to be capable to operate safely within 
their flight envelope, at higher angles of attack, at superior speeds, and without weight penalties. 
The implementation of smart material systems could play an important role in the design of 
future advanced space vehicles. 

As is well known, in the determination of both the dynamic response to time-dependent 
excitations and of resonant conditions as well as in the flutter analyses, the natural frequencies are 
an important physical parameter that intervenes in an explicit way [R2]. For this reason, one of 
the features of adaptive structural technology applied to aeronautical structures is the ability to 
conveniently control the eigenfrequencies ofthe system. In an effort to contribute to this 
problem, the free vibration analysis of adaptive aircraft wing structures incorporating induced 
strain actuation [R3,4] will be considered in this paper. The wing structure modelled as a 
thin— walled beam is assumed to be composed of an induced strain actuator, e.g. a piezoelectric 
layer superposed on the master structure. 

The global constitutive equations associated with a piezoelectric material layer (the 
induced strain actuator), polarized in the thickness direction and exhibiting transversely isotropic 
properties, superposed on the thin-walled beam structure also made of a transversely-isotropic 
material, are derived. These equations, used in conjunction with a generalized Hamilton’s 
variational principle [R5,6], yield the equations governing the motion of cantilevered adaptive 
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wing structure modelled as a thin-walled beam. Coupling the basic properties of piezoelectric 
materials (i.e., their actuating capabilities) with a control law, it will be shown that the natural 
frequencies can be controlled in a known and predictable manner. 


BASIC ASSUMPTIONS 

The thin— walled beam model used in the present paper is based on the following kinematic 
statements [R8— 12]: 

1) The original shape of the cross-actions of the beam is preserved. 

2) The transverse shear flexibility exhibited by the advanced composite material 
structures is in co rp orat e d . 

3) The non-uniform torsional model is considered. In this respect, the rate of twist 
assumed to vary along the beam axis constitutes a measure of the waxping restraint 
effect. The primary warping displacement, throughout the cross-section of a beam, 
is assumed to have a similar distribution as the one associated with the St. Venant 
(uniform) torsion theory. 

4) Incorporation of the secondary warping [R9] whose effect for composite material 
structures ccmld be comparable to the primary warping. 

5) In addition to these statements of kinematic nature, another one, of a static nature 
is adopted [RIO]. According to this statement, the hoop stress resultant N gB is 

considered negligibly small with respect to the remaining ones. 

Based on the above assumptions, the displacement field can be expressed as [see R7-9]: 



u( x »y»®) — u Q (z) — y0(z) , 

(1) 


v(x,y,z) = v Q (z) + x0(z) , 

(2) 


w(n,s,z) - w 0 (z) + x(s)# y (z) + y(s)* x (z) -FJs)©'(z) 



+ n [a£ V') • 

( 3 ) 

where 



«,(■) = 

( 4 ) 


» x (*) = V (, ) -Y o w ’ 

( 5 ) 



(6) 

The warping function is expressed as 



Fj8)*/ 8 [r n (8)“lflds, 

( 7 ) 
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( 8 ) 


where the torsional function j> is given by 


and 



( 9 ) 


As a result, six kinematic variables Uq(z), Vq(z), Wq(z), 0 (z), 0 x (z), and 9(z) representing three 

translations in the x, y, z directions and three rotations about the y, x, and z axis, respectively, 
are used to define the displacement vector (i.e., the displacement components u, v, and w in the x, 
y, and z directions, respectively). Here (s,z,n) and (x,y,z) will be used to denote the surface and 
cross-section reference coordinate systems, respectively (see FI). The quantity h[= h(s)] denotes 
the wall thickness of the beam (allowed to vary along the periphery); Aq denotes the 

cross— sectional area bounded by the mid-line; 0 denotes the total length of the contour mid-line; 

| (• )ds denotes the integral around the entire periphery C of the mid-line cross-section of the 
J C 

s 

beam; while / r (s)ds[=fi(s)] is referred to as the sectorial area, 
o 

Based on the kinematic representations, Eqs. (1)— (6), the strain measures assume the form: 
Axial Strain: 


S zz ( n ,s, z ) = 5 zz (s,z) + n S zz (s,z) , (10) 

where 

5 zz (s,z) = w'(z) + 0y(z)x(s) + 0'(z)y(s) - ©' \z)Fjs) , 
and (11) 

Kz( b >*) = gf “ d x ( z ) af~ e ( z ) a ( 8 ) > 

are the axial strains associated with the primary and secondary warping, respectively. 

Membrane Shear Strain: 



Transverse Shear Strain: 

s m (m) = [ 0 y M + «o(*)l aJ - t« x W + »oWl a? ( l3 ) 

Within the present theory the warping measure is expressible as 
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w M = e'(i) . 

Here, and in the following developments (•) =d(')j dz. 


( 14 ) 


PIEZOELECTRIC CONSTITUTIVE EQUATIONS 


The linear constitutive equations for a 3— D piezoelectric continuum expressed in the 
contracted indidal notations are (see e.g., R13,14): 


a '\ ~ C ij S j e ld 3c’ 
D k = e kj S j + h ’ 


(15) 


where <r. and Sj (i j = T^5) denote the stress and strain components, respectively, where 



Spr when p = r , j = 1, 2, 3 
2S pr when p * r , j = 4, 5, 8 


C jj, 6j^, are the das tic (measured for conditions of constant electric field), piezoelectric, and 

dielectric constants (measured under constant strain), while and (k = 1^3) denote the 

electric field intensity and electric displacement vector, respectively. In Eqs. (15) the summation 
over repeated indices is implied. While Eq. (15)j describes the converse piezoelectric effect 

( consisting of the generation of mechanical stress or strain when an electric field is applied), Eq. 
( 15)2 describes the direct piezoelectric effect (consisting of generation of an electrical charge when 

a mechanical force or pressure is applied). 


In a piezoelectric adaptive structure the direct effect is used for distributed sensing while 
the converse effect is used for the active distributed control. Equations (15) are valid for the most 
general anisotropic case, i.e., for tridinic crystals. In the following, we will restrict the generality 
to the case of a transversely isotropic continuum, the n-axis being an axis of rotary symmetry 
coinciding with the direction of polarization (thickness polarization). 

In this case the piezoelectric continuum is characterized by 5 independent elastic constants, 
3 independent piezoelectric constants and 2 independent dielectric constants. In the following 
developments, for the sake of simplicity, we will assume that the master structure is constructed 
from transversely isotropic material layers, the axis of symmetry bring parallel to the n-axis. 

We will also assume that the riectric field vector ^ is represented in terms of its 
component ^ only, (implying ^ = 0) and, as a result of the uniform voltage distribution, 

is independent of space (but, possibly on time). In matrix form the constitutive equations are: 




where the index k in brackets affecting a quantity identifies its affiliation to the k-th layer. 
In terms of the engineering constants the coefficients C- • associated with a 

transversely— isotropic continuum are expressed as (where for the sake of simplification the 
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superscript 9 was discarded): 

C u = (Ei/ 2 -e')E/A, 

'j 

C 12 = -(E 1 /' 2 + e' I /)E/A, 

— — u (l+i/)EE /A , (18) 

C m = -(1-^)e' 2 /A , = G , 

^1 1*^12 p/_ E \ 

modulus and 
, respectively while 

G' denotes the transverse shear modulus. Equation (16) reveals that a piezoelectric continuum 
exhibiting anisotropic behavior cannot generate or detect shear stresses/strains by applying or 
detecting electric fields along the n-direction. 

INCORPORATION OF ACTUATOR PATCHES AND ASSOCIATED 
CONSTITUTIVE EQUATIONS 

As was stated earlier, the master structure is composed of layers of elastic materials also 
exhibiting transversely— isotropic behavior. Their associated constitutive equations could be 
established formally by discarding the electrical effects in Eq. (16) and Eq. (17). 

Suppose that the master structure is composed of such layers, while the actuator 
(superposed on the master structure) is composed of l piezoelectric layers. We also stipulate that 
the actuators are distributed over the entire span of the wing (i.e., along the entire span wise 
coordinate z), (see F2) while along the circumferential s- and transversal n-directions they are 
distributed according to the law (see F3): 

R( k )(n) = H(n - n ) - H(n - 

R( k )(») = H(s - s ) - H(s - 

where H denotes Heaviside’s distribution and R is a spatial function [R15]. In terms of the 
coordinates (s,z,n), related to the beam, the constitutive equations are expressed as 


n . ) , 
(k + ) 

8 , ) , 
(k + ) 


(19) 


where A = (1 +p) (2E*/ 2 + E v - E ) while E, i/, and E . v denote Young’s 
Poisson’s ratio in tne plane of isotropy and transverse to tne plane of isotropy 
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4? 4 k) R (k)(“) R (k)W 


(k) 

(k) = c ^k) s (k ) 


4? 4 k) R (k)( n ) R (k)(®) 


nz 


nz 


( 20 ) 


LOCAL BEAM CONSTITUTIVE RELATIONSHIPS 

Assume that the master structure is composed of m layers while the actuator part is 
composed of l layers. As a result, the global beam stress resultants and couples could be obtained 
through the integration of the 3-D stress components through the laminate thickness and 
afterwards through the integration of the local stress resultants along the contour of the beam. 
Invoking the assumption 5), the local beam stress resultants defined in terms of the associated 
strains are: 


Stress Resultants : 


N z* = R 11 s zz + R n S Z, 
N — A/./. S < 

8Z 00 SZ ’ 

N = A aa S , 
nz 44 nz 

Stress Couples: 

I'„ = ^n S w + K 11 § 
zz IX zz 11 zz 

sz 66 sz 


_ N a 
11 ** > 


T a 

> 


( 21 ) 


In these equations K.. denote the modified local stiffness coefficients while N* z and L zz the 
piezoelectrically induced stress resultant and stress couple. Their expressions are 


R ll~ A 11 


12 


) 


'll * B 11 


B 12 A 12 
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(22) 


whereas 





B 

X 


2 

12 

11 




l 

L «= l 4 k) 4i )R (k) ( *'*)k.+, 

k=l ** > 




(23) 




m +1 

n c i ( n (k) ~ tt ?k-i)] 

k=l 


m+f 


-n cjs>h 


k=l 


3 

(k) 



(24) 


define the stretching, bending-stretching and bending stiffness quantities, respectively. In Eqs. 
(23), for the sake of brevity the following notation 

R (k)( 8 » z ) 5 ^(k)^(k)( 1 ^ (^®) 

was introduced. 

It should be emphasized that the constitutive equations (21) relate the stress resultants and 
couples with the strain components and the electric field, The analysis yielding such 

equations is called uncoupled because an approximate solution to the electrostatic problem is 
postulated, with the result that is constant. 


THE DYNAMIC EQUATIONS OF ADAPTIVE THIN-WALLED BEAMS 


In the previous sections the kinematic equations and the local constitutive equations for 
the considered adaptive structure were derived. In order to obtain the equation of motion of an 
adaptive TWB and the associated boundary conditions, we will make use of Hamilton’s 
variational principle extended to the case of a linear 3-D piezoelectric medium [R5,6]: 



dt [[ [jp UjUj - <#+ f.U.] dr + J n (* k flJ k - 


a6<t>) dfl] = 0 


(26) 
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where ^denotes the electric enthalpy density defined as: 




( 27 ) 


r and fl denote the volume and bounding surface of the continuum, respectively; denotes the 

specified surface-fraction; a the applied surface charge; ij the body forces; tg, t j denote two 

arbitrary instants of time; p the material density; 6 the variation sign, while the overdots denote 
time derivatives. 

Recalling that 

|^=a. and = “ D k ’ ( 28 ) 

when = 0, defining in terms of a potential function as 

<$ = -*i . W 

taking the variations in (26) and applying Green’s theorem whenever possible, we obtain, by 
considering the variations dTJ., 6<f> and tfa. as independent and arbitrary, a known version of the 

fundamental equations of the 3— D linear piezoelectricity theory, namely 





in t , 


ff i = 


»J J 


<f> = $ on ft 


(30) 


Uj = U- a- njD. on fl u , 

are obtained, where n. denote the components of the normal to 0. Consideration in Eq. (26) of 

the displacement components U- as u Q , v Q , w Q , $ x , 9^, and 9, and of Sj as the ones defined by 

Eqs. (10) — (13); employment of dr = dndsdz and performing the integration across the s— and 
n— directions, the 3— D problem is reduced to an equivalent 1— D one, in which all the quantities are 
dependent on the z— coordinate only. To this end, the 1— D stress resultants and stress couples 
specific to the theory of TWBs defined in terms of the local (i.e., 2— D) beam stress— resultants are 

T ,« 5 f c V; <y«) = f c K a? + N *> &] d » • 

=f c K^ _N “^ d ’ • M x (,) - - 

(31) 
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and 


2A 

M yW = | c (^zz +L I .S] d * :M » W “f c N »ir d,t 
B JC-) = $ c PJfiK m + »(‘)I. Jd. . 

These quantities are wfeticd to as the axial and shear forces (in the x and y directions), and the 
mome n ts (in the x, y and s directions) as well as tbs bimoment global quantities. In trams of the 
1— D stress resultant measures, (Eq. 31), the equations of motion of adaptive TWBs become: 


*>o : Q x _I l + p x = 0 > 

*v VV-p y =°. 

T » -I 3 + Ps = 0 ’ 

B L' +M z-( I 4- I 8 / ) + m z + b w= 0 - 
M x : M x-Q y -I 5 + m x = ° > 

S9 f M y“Qx“ I 7 + m y = 0 * 


(32) 



•-directions, respectively; b^ is the bimoment of the external loads while Ij (i = I^ff) are the 

inertia terms not displayed in the paper. By virtue of Eqs. (21), (23) and (31) the stress resultant 
T z , the stress couples M x and M y as well as the bimoment could be recast in a form in which 

the actuator effect appears in a separated form, namely 



M x” V*x- 


B =B -B , 
fal u u * 


(33) 


where the quantities affected by an overhat (*) and an overtilde (~) identify the pure mechanical 
and piezoelectric contributions to the indicated quantities, respectively. 

Being concerned in this study with the free vibration problem only, the loading terms 
occurring in Eqs. (32) may be discarded. From the Hamilton’s principle, in addition to the 
equations of motion (32), the boundary conditions (BCs) are also obtained. For the case of the 
beam clamped at z = 0 and free at z = L, the BCs are: 
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M z + B >¥z; B u =? w - <*> 

where the undertilde sign affects the prescribed quantities and where L denotes the length of the 
beam. It could be verified that consistent with seven boundary conditions at each edge, a 
fourteenth order governing equation system is obtained. 

In the case of the general anisotropy of the layer materials (i.e., of the master structure, of 
the actuator patches or of both of them), the system of governing equations results in a complete 
coupled form. However, in the present case of anisotropy, the governing equations obtainable 
from Eqs. (32)^ g appear decoupled from the ones obtainable from (32)2 5 and (32)g ^. In other 

words, in the present case of anisotropy, the governing system of equations and the associated BCs 
are splitting exactly into three uncoupled sub— systems. While the first mentioned sub— system of 
equations governs the lateral bending motion (or in helicopter terminology, the flap— lag motion), 
the second and the third ones govern the vertical bending and twist motions, respectively. 

In this study only the undamped free bending vibration case is studied. The pertinent 
governing equations for free vibration of thin— walled beams are: 

For transverse bending vibrations: 


*55^’ + O + vS. = 0 

*33*x “ a S5 ( T o + B x) + “^ b 4 + b 14^x = °' 

(36) 

with the BCs at z = 0: 


0 

II 

II 

>° 

(37) 

and at z = L: 


V +e x = 0 and a 33*x = *x 

(36) 

For lateral bending vibrations: 


+ G y) + h l^ n o = 0 

(39) 
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( 40 ) 


*22 V ’ " a 44(“o’ + V + = 0 

with the BCs at z = 0: 



and at s — L: 


\ +*y = 0 ** d *22 V = ®y ( 41 ) 

These equations were obtained by considering the following harmonic time-dependence 

^ r (s,t)= 3 r (z)exp(iut) , i = (42) 

where ^stands for a generic field variable. 

Both the governing equations and the associated BCs art expressed in terms of the 
quantities involving the dependence on the »-coordinate only, i.e., according to Eq. (42) in terms 
of the barred counterparts of the field variables. However, for the sake of simplicity these 
overbars have been dropped. The expressions of elastic constants a^ and mass terms b. 

intervening in Eqs. (36) and (39) are not displayed here. In the same equations M x and M y 
denote the piezoelectrically induced moments expressed as: 

*x = i c J i 4 k) [- (k+) - l 1 -sit) 



For the case of the symmetrically located actuators across the thickness of the beam structure, the 
terms effected by a solid line in Eqs. (43) vanish. In light of the actuator configuration, it may be 
inferred that the expressions (43) are independent on the z—coordinate. This explains why their 
contribution in the governing equations is immaterial and why in the BCs they intervene as 
nonhomogeneous terms. 
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THE CONTROL LAW 

The adaptive nature of the wing (or rotor blade) structure is introduced by requiring the 
applied electric field ^ to be dependent to one of the mechanical quantities of the structure in 

motion. With this in mind, a number of control laws could be implemented. Their efficiency may 
be measured by their ability to modify as much as possible the natural frequencies of the structure 
with respect to the energy input. In the present considerations, two control laws were considered. 

The two independent control laws require that: i) the applied electric field ^ is 

proportional to the vertical or lateral bending moments at the wing root, M x (o) or M^(o), 

depending on whether the control of the natural frequencies is associated with the vertical or 
lateral motion, respectively and ii) the applied electric field ^ is proportional to the vertical 

( V 0 (L)) or lateral (u Q (L)) deflections of the beam tip, depending on the two directions of 

vibrations whose frequencies are to be controlled. In light of these control approaches, we may 
formulate the following control laws (labelled as CL1 and CL2) by involving also the boundary 
conditions given in expression (38) and (41). 

CL1) 

a) Associated with the vertical bending 

» x '(t) = V x '(o) («) 

b) Associated with the lateral bending 

«/(L) = V,'(o) ( 45 ) 

CL2) 

a) Associated with the vertical bending 

d x '(L) = R p v 0 (L) (46) 

b) Associated with the lateral bending 

«y'(L) = 8 p u 0 (L) (47) 


As can be seen from Eqs. (44)-(47), the feedback gains and 


K p are dimensional and 


nondimensional, respectively. The nondimensional counterpart 

K_ = KL 2 
P P 


of R is R defined as 
P P 
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These control laws express the fact that the control bending moments at the wing tip 
(induced by the piezoelectric strain) should be proportional as per the CL1 to the mechanical 
naniBng moments at the wing root, and within the CL2 to the transversal or lateral deflections at 
the wing tip. 


NUMERICAL ILLUSTRATIONS AND DISCUSSION 

The adaptive wing is modelled as a symmetric composite box— beam. In order to control 
the vertical bending frequencies, the piezoceramic actuator layers (selected to be of PZT— 4 
ceramic) are located on the upper and bottom surfaces of the master structure only, while in order 
to control the lateral bending frequencies the piezoceramic layers are located on the lateral walls 
of the master structure. The geometrical characteristics (restricted to the former situation, only) 
are displayed in F2 while the constants characterizing the PZT— 4 piezoceramic (see [R16]) are: 


C n = 2.016 X 10 7 psi; 

C 12 = 1.128 x 10 7 psi 

C 13 = 1.0776 x 10 7 psi; 

C 33 = 1.6679 x 10 7 psi 

C u = 3.7128 x 10 6 psi; 


e^ = -0.0297 lb/in V; 

e 33 = 0.08623 lb/inV 


e lg = 0.07252 lb/inV. 

For the of simplicity assume that the master structure is composed of a 
transversely— isotropic material whose elastic characteristics are identical to the piezoceramic 
actuator layers. 

By using the exact approach devised in [R17— 19] and extended afterwards in [R7-9], the 
eigenvalue problems associated with the vertical and lateral vibrations of TWBs yield the 

V V 

variation of the nondimensional eigenfrequendes u vs. the feedback gains and K^. 

In F4-7, by using the control law CL1, the variation of the first four eigenfrequendes 
assodated with the vertical bending vs. the feedback gain was diagrammaticauy represented, 
while in F8 and 9, by using the control law CL2, the variation of the vertical (plunging) and 

lateral bending fundamental frequences u> vs. the feedback gain was obtained. F4— 7 reveal 

that within control law CL1 the increase of the odd eigenfrequendes occurs for negative feedback 
gains while the increase of the even eigenfreqnendes requires implementation of positive feedback 
gains. F8 and 9 (as well as the results obtained, but not displayed in the paper) reveal that within 
control law CL2 the increase of the eigenfrequendes occurs generally for positive feedback gains. 
Within a CL2 type control a similar trend was also obtained in R20. 

It shoula be pointed out that in all diagrams, F4-9, the frequencies were normalized with 
respect to the ones corresponding to the non— adaptive structure (for which case K = 0). 


Having in view that, roughly speaking, w ~ D*/ ~ h^, where D and h are the associated 
bending rigidity and thickness of the beam, we could infer that the linear increase of 
eigenfrequenties would have been accomplished without thehdp of this adaptive technology 
through an unaffordable weight penalty. The obtained results reveal, once again, the great 
importance of the implementation of adaptive technology applied to the control of the material 
frequencies of the structure. In short, the obtained results reveal that by using the adaptive 
properties of the structure it is possible to increase the eigenfrequendes of the system and 
consequently to modify in a beneficial way the dynamic response characteristics of the structure, 
whose role in the present case could be played by an aircraft wing or a rotor blade structure. 
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Fig. 3 Actuator patch. 
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Fig . 4 Variation of the normalised plunging eigenfrequency vs. the control gain (first 

mode and first control law). 
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Variation of the normalized plunging >'igenfrequency vs. the control gain 
(second mode and first control law). 
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Fig. 6 Variation of the normalized plunging eigenfrequency vs. the control gain (third 
mode and first control law). 
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Variation of the normalised plunging eigenfrequency vs. the control gain (first mode 
and second control law). 
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